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We study ferrimagnetism in the ground state of the antiferromagnetic Heisenberg model 
on the spatially anisotropic kagome lattice, in which ferrimagnetism of the conventional Lieb- 
Mattis type appears in the region of weak frustration whereas the ground state is nonmagnetic 
in the isotropic case. Numerical diagonalizations of small finite-size clusters are carried out to 
examine the spontaneous magnetization. We find that the spontaneous magnetization changes 
continuously in the intermediate region between conventional ferrimagnetism and the nonmag- 
netic phase. Local magnetization of the intermediate state shows strong dependence on the site 
position, which suggests non-Lieb-Mattis ferrimagnetism. 
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Ferrimagnetism has been studied extensively as an im- 
portant phenomenon that has both ferromagnetic nature 
and antiferromagnetic nature at the same time. One of 
the fundamental keys to understanding ferrimagnetism 
is the Marshall-Lieb-Mattis (MLM) theorem. 1 ' 2 ) This 
theorem clarifies some of the magnetic properties in the 
ground state of a system when the system has a bipar- 
tite lattice structure and when a spin on one sublattice 
interacts antifcrromagnetically with a spin on the other 
sublattice. Under the condition that the sum of the spin 
amplitudes of spins in each sublattice is different between 
the two sublattices, one finds that the ground state of 
such a system exhibits ferrimagnetism. In this ferrimag- 
netic ground state, spontaneous magnetization is realized 
and its magnitude is a simple fraction of the saturated 
magnetization. We hereafter call ferrimagnetism of this 
type the Lieb-Mattis (LM) type. 

Some studies in recent years, on the other hand, re- 
ported cases when the magnitude of the spontaneous 
magnetization of the ferrimagnetism is not a simple frac- 
tion of the saturated magnetization. 3 ~ 9 ) The ferrimag- 
netic ground state of this type is a nontrivial quantum 
state whose behavior is difficult to explain well only 
within the classical picture. Ferrimagnetism of this type 
was first predicted in ref. 10 using the quantum rotor 
model. The mechanism of this ferrimagnetism has not 
been understood sufficiently up to now. Hereafter, we call 
this case the non-Lieb-Mattis (NLM) type. Note that in 
the cases when NLM ferrimagnetism is present, the struc- 
ture of the lattices is limited to being one-dimensional. 
Recall that the above conditions of the MLM theorem 
do not include the spatial dimension of the system; the 
MLM theorem holds irrespective of the spatial dimen- 
sionality. We are then faced with a question: can NLM 
ferrimagnetism be realized when the spatial dimension is 
more than one? 

The purpose of this letter is to answer the above ques- 
tion concerning the existence of NLM ferrimagnetism 



in higher dimensions. In this letter, we consider a case 
when we introduce a frustrating interaction into a two- 
dimensional lattice whose interactions satisfy the con- 
ditions of the MLM theorem. When the frustrating in- 
teraction is small, ferrimagnetism of the LM type sur- 
vives; however, the ferrimagnetism is destroyed with the 
increase in the frustrating interaction and the system 
finally becomes nonmagnetic due to the considerably 
large frustrating interaction. We examine the behavior 
of the collapse of the ferrimagnetism and the existence 
of an intermediate region between the LM ferrimagnetic 
and nonmagnetic phases by means of the numerical- 
diagonalization method applied to finite-size clusters. 
Our study of the two-dimensional system successfully 
clarifies the existence of the intermediate phase and cap- 
tures a feature of NLM ferrimagnetism. 

First, we explain the model Hamiltonian examined in 
this letter. The Hamiltonian is given by 
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where Si denotes an S = 1/2 spin operator at site i. 
Sublattices A, B, and B' and the network of antiferro- 
magnetic interactions J\ and Ji are depicted in Fig. 1. 
Here, we consider the case of isotropic interactions. The 
system size is denoted by N s ; the saturation magneti- 
zation is M sat — N s /2. Energies are measured in units 
of Jx\ thus, we take Ji = 1 hereafter. We examine the 
properties of this model in the range of < J2/J1 < 1. 
Note that in the case of J2 = 0, sublattices B and B' 
are combined into a single sublattice; the system satis- 
fies the above conditions of the MLM theorem. Thus, 
ferrimagnetism of the LM type is exactly realized in this 
case. In the case of J2 = Ji, on the other hand, the 
lattice of the system is reduced to the kagome lattice. 
The ground state of the system on the kagome lattice 
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Fig. 1. Network of interactions in the system and sublattices A, 
B, and B'. Black straight lines and green dotted lines denote 
interactions of Ji and J2, respectively. Open circles at lattice 
points represent 5=1/2 spins. The system of classical spins in 
this lattice was studied by Monte Carlo simulations in ref. 11. 

without a magnetic field is known to be singlet from 
numerical-diagonalization studies, 12 ~ 15 ) which indicates 
that the ground state is nonmagnetic. One thus finds 
that LM ferrimagnetism collapses between J2 = and 
J2 = J\. Consequently, we survey the region between 
the two cases. 

Next, we discuss the method we use here, which is nu- 
merical diagonalization based on the Lanczos algorithm. 
It is known that this method is nonbiased beyond any ap- 
proximations and reliable for many-body problems such 
as the present model. A disadvantage of this method is 
that the available system sizes are limited to being small 
because the dimension of the matrix grows exponentially 
with respect to the system size. To treat systems that are 
as large as possible, we have developed parallelization in 
our numerical calculations using the OpenMP and MPI 
techniques, either separately or in a hybrid way. 16 -* 

In this letter, we treat the finite-size clusters depicted 
in Fig. 2 when the system sizes are N s = 12, N s = 24, 
N s = 27, and N s — 30 under the periodic boundary con- 
dition and iV s = 33 under the open boundary condition. 
Note that each of these clusters forms a regular square 
although clusters (b) and (d) are tilted. The next larger 
size under the condition that a regular square is formed 
is N s = 48, which is too large to handle using the present 
method, even when one uses modern supercomputers. 

Before our numerical-diagonalization results for the 
finite-size clusters are presented, let us consider the direc- 
tions of the spins in the ground state within the classical 
picture. We here examine the spin directions of classi- 
cal vectors with length S depicted in Fig. 3. One ob- 
tains the energy of the spin state with angle 9 to be 
E/J x = -(2N s /3)S 2 [2cos(ir-9) + {J 2 /J 1 )cos(26)}. This 
expression of the energy indicates that for J2/J1 < 1/2, 
the state of 9 = 0, namely, ferrimagnetism of the LM 
type, is realized. Thus, the normalized magnetization of 
this state is M/M sa t = 1/3. One finds, on the other 
hand, that for J2/J1 > 1/2, the lowest-energy state is 
realized for nonzero 9 when J1/J2 = 2cos(#) is sat- 




Fig. 2. Finite-size clusters: (a) N s = 12, (b) Ns = 24, (c) N s = 
27, and (d) N s = 30 under the periodic boundary condition, 
where red dashed lines denote a single finite-size cluster with each 
system size. Black straight lines and green dotted lines are the 
same as in Fig. 1. Note that cluster (c) under the open boundary 
condition includes N B = 33 spins. 
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Fig. 3. Ferrimagnetic spin direction in the classical picture. 

isfied. The normalized magnetization of this state is 
M/M sat = - l)/3. When J 2 /Ji becomes unity, the 
magnetization finally vanishes. This classical argument 
will be compared with our finite-size results obtained 
from numerical diagonalizations. 

Now, let us present our numerical results for the quan- 
tum case. First, we show our data for the lowest energy in 
each subspace of 5*'°*, which reveal the magnetization of 
the systems. Figure 4 depicts our results for the system 
with N s = 30 depicted in Fig. 2(d). Note that M sat = 15 
in this case. For J2/J1 = 0.5, the energies from S"* ' = 
to 5*°' = 5 are numerically identical, which means that 
M/M sa t becomes 1/3 and that ferrimagnetism of the LM 
type is realized. For J2/J1 = 1, the energy for Sl ot = is 
lower than the other energies for larger S* ot . The ground 
state of this case is nonmagnetic. For J2/J1 = 0.6, the 
energies from Sl ot = to Sl ot = 2 are the same; thus, 
we find that the spontaneous magnetization is M = 2, 
which is smaller than the value for ferrimagnetism of the 
LM type. One finds that a state with intermediate mag- 
netization appears between LM-type ferrimagnetism and 
the nonmagnetic state, at least according to the finite- 
size calculations. 

Next, we examine the region of such an intermediate 
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Fig. 4. Lowest energy in each subspace of S* ot for the system of 
7V S = 30 depicted in Fig. 2(d). Results for J2/J1 = 1, 0.5, and 
0.6 are presented by black circles, blue triangles, and red squares, 
respectively. Inset: our data in the entire range of S!; ot given for 
Ja/Ji = 1 and 0.5. 



state for various system sizes; our results are depicted in 
Fig. 5. In the case of N s = 12, the intermediate state 
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condition is close to that for the N s = 24 case under the 
periodic boundary condition. This is consistent with the 
fact that there are 21 sites in the inner part of cluster 
(c) of N s = 33 under the open boundary condition. Our 
present results for both boundary conditions imply that 
the presence of the intermediate phase is irrespective of 
the boundary conditions. 

An important characteristic of NLM ferrimagnetism is 
that the local magnetization in sublattice exhibits long- 
distance periodicity, which is absent in LM-type ferri- 
magnetism. Note that one cannot detect this periodicity 
in the cases under the periodic boundary condition. We 
thus examine the local magnetization in the intermedi- 
ate phase for the case under the open boundary con- 
dition; the results for N s = 33 are depicted in Fig. 6. 
For J2/J1 = 0.5 with LM-type ferrimagnetism, the local 
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Fig. 6. Local magnetizations of the cluster with N s = 33 under 
the open boundary condition. Results for J2/J1 = 0.5, 0.53, and 
0.57 are presented for the sites surrounded by the blue rectangles 
in the inset. Site numbers 1 to 6 correspond to the sites from left 
to right in the blue rectangle. 



Fig. 5. Dependence of the spontaneous magnetization normalized 
by the saturated magnetization on J1/J2. The solid line repre- 
sents the result for the magnetization within the classical picture 
shown in Fig. 3. Note that we take the J1/J2 dependence as the 
abscissa because the classical magnetization shows linear depen- 
dence not on J2/J1 but on J1/J2. Results for N s = 12, 24, 27, 
and 30 under the periodic boundary condition are presented by 
black pluses, violet crosses, green squares, and blue diamonds, 
respectively. Red circles denote results for N B = 33 under the 
open boundary condition. 

between LM-type ferrimagnetism and the nonmagnetic 
state is absent; on the other hand, the intermediate re- 
gion exists for all the larger systems. Note that the width 
of the intermediate region increases for the cases under 
the periodic boundary condition when N s is increased. 
This indicates that the intermediate phase is present in 
the thermodynamic limit. One of the characteristics ob- 
served is that the continuity of the magnetization im- 
proves with increasing N s . In the cases under the open 
boundary condition, we successfully detect the interme- 
diate phase although its width is relatively smaller. The 
width for the N s = 33 case under the open boundary 



magnetization shows weak dependence on the position 
of sites, although (Sf) at edge sites 1 and 6 is slightly 
larger than those at interior sites, where the site numbers 
are illustrated in the inset of Fig. 6. This small difference 
originates from the edge effect due to the open bound- 
ary condition. For J2/J1 = 0.5, the edge effect does not 
seem to affect (Sf) at internal sites. For J%/J\ — 0.53 
and 0.57, on the other hand, (Sf) at edge sites 1 and 6 
becomes very small. The behavior of this appearance of 
the edge effect is different from the case of JijJ\ = 0.5. 
For J2/J1 — 0.53 and 0.57, one finds a strong dependence 
of (Sf } on the position of the site from site 2 to site 5. For 
J2/J1 = 0.53, (Sf) at sites next to the edges seems to be 
affected by the edge sites. It is unclear whether or not the 
case of J2/J1 = 0.53 corresponds to NLM type ferrimag- 
netism at present. For J2/J1 = 0.57, on the other hand, 
the strong dependence on the site position suggests the 
existence of origins that are different from the edge effect. 
The system size N s = 33 is sufficiently small for long- 
distance periodicity to be observed clearly. Although the 
present results are not decisive evidence of the periodic- 
ity, our finding of the large change in (Sf) is considered 
as possible evidence. In order to obtain decisive evidence, 
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calculations on systems of larger sizes are required, which 
are unfortunately difficult at the present time. Instead of 
the present two-dimensional kagome case, we are now 
examining a quasi-one-dimensional system on a kagome 
stripe lattice. Both systems partly share the same lat- 
tice structure. The system on the kagome stripe lattice 
reveals the clear appearance of NLM ferrimagnetism in 
the intermediate region. 18 ) Results will be published else- 
where. 

The phenomenon of ground-state magnetization 
changing continuously with respect to a continuous pa- 
rameter in a model Hamiltonian has been reported in 
other cases. Tonegawa and co-workers reported such a 
phenomenon in spin systems with anisotropic interac- 
tions. 19 ~ 23 ) It is unclear at present whether or not the 
states of this continuously changing magnetization in the 
anisotropic case show long-distance periodicity because 
the behavior of local magnetization has not been inves- 
tigated yet. Since this phenomenon disappears in the 
isotropic case when the quantum effect is stronger than 
that in the anisotropic case, this phenomenon is con- 
sidered to arise from the anisotropy. From this point of 
view, the origin of this phenomenon seems to be different 
from that of intermediate ferrimagnetism in the isotropic 
case studied here. Another reported phenomenon is par- 
tial ferromagnetism in the Hubbard model 24 ' 25 ) when the 
system is hole-doped near the half-filled Mott insulator. 
The origin of this phenomenon has been clarified to be 
the formation of spin polarons around doped holes. The 
mechanism of these two cases is different from the present 
case of NLM ferrimagnetism. 

Finally, we briefly discuss possible future experiments. 
For volborthite, eq. (1) was proposed as a model Hamil- 
tonian from the argument of its crystal structure, 26 ' 27 ) 
although NLM ferrimagnetism has not yet been ob- 
served in this material. A theoretical study on the spatial 
anisotropy of this material indicated that the deviation 
of the anisotropy from the isotropic kagome point is not 
particularly large. 28 ) This is consistent with our present 
result because the nonmagnetic ground state is realized 
around the region of weak anisotropy as shown in Fig. 5. 
In order to observe NLM ferrimagnetism experimentally, 
it is necessary to realize a case with larger anisotropy. 
The measurement of volborthite under high pressure in 
the direction of the a-axis or discoveries of new materials 
might lead to such an observation. 

In summary, we have clearly shown the existence of a 
ground state of non-Lieb-Mattis type ferrimagnetism in a 
two-dimensional lattice that lies between the well-known 
Lieb-Mattis type ferrimagnetic phase and the nonmag- 
netic phase including the kagome-lattice system. The 
nontrivial ferrimagnetism we have found in the interme- 
diate phase occurs as a consequence of magnetic frustra- 
tion. Our present result indicates that non-Lieb-Mattis 
ferrimagnetism is a general phenomenon irrespective of 
the spatial dimensionality. 
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